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I^IBGEDING NOT FH.MBD 


1.0 SUMMARY 


The purpose of this study was to review data and theory pertinent 
to the autoignition of liquid oxygen/liquid hydrogen {LOX/LH 2 ) and 
liquid oxygen/RP-1 (LOX/RP-1) propellants and develop, where possible, 
physical models of the processes supporting or contributing to auto- 
ignition. The pertinent data included that developed in a study led 
by Dr. Erich Farber (NASA Contract NAS 10-1255)^ and other preceding 
explosive tests such as Project PYRO.^ During the Farber study, a 
concept called the "Critical Mass Hypothesis" was developed which 
stated that when fuel and oxidizer are mixed with a certain energy, 
there is a mixed mass where explosion is 100% probable, called the 
"Critical Mass". The hypothesis is based on two simultaneous pro- 
cesses: 1) mixing of the propellants, and 2) development of electro- 

static charge sufficient to cause ignition sparks. This is known as 
"autoignition". Farber 's work was aimed solely at getting empirical 
techniques for predicting blast yield (TNT equivalents) for liquid rocket 
propellants. However, the objective of this study was to describe the 
physical environment as fully as possible and relate this description 
to the autoignition phenomenon. 

Processes other than electrostatic charge were also explored as 
possible alternative or competing ignition mechanisms. Two in particular 
were cosmic rays and heat from ortho-parahydrogen conversion. It was 
found that the energy levels developed from cosmic ray ignition, while 
measurable, is far less than that needed to ignite an explosive 
hydrogen-oxygen mixture. Sufficient total energy is available in the 
ortho-parahydrogen conversion (from 95 to 99% para) to cause ignition, 
but the temperature is likely to be too low and the heat carried off 
easily by the surrounding fluid. Therefore, neither of these processes 
is considered a significant ignition mode when compared to electrostatic 
phenomena. External ignition sources (such as sparks from tearing metal) 
are not necessarily present and as such have low probability as compared 
with internal processes associated with autoignition by static electrical 
discharge. 

A mixing model was developed which gave a relationship between mass 
which is mixed and mixing energy for energies above the boiling energy. 
The relationship, developed from first principles, verified the form of 
the high energy term in the critical mass equation. 

The study of static electrical effects included three main topics. 
The first involved estimating the streaming potentials built up between 
the liquids during mixing. It was demonstrated that a sufficient elec- 
trical field for ignition, as predicted by streaming potential theory, 
could be built up by the relative motion of the fluids. The second was 


concerned with discharge within the liquid. The dielectric breakdown 
within the liquid appears to be a plausible mechanism for ignition since 
discharge across a 0 . 06 -cm space through a stoichiometric mixture would 
be a reasonable ignition source. The third considered build-up and 
discharge at the surface and the speculation concerning probability of 
ignition at the surface of the liquid mixture. It is plausible that 
droplet charging and droplet-droplet discharge at the surface could cause 
autoignition. In all these topics specific numbers could not be devel- 
oped, but parameter range applicable to the cases studied were of the 
order of magnitude such that electrostatic discharge was reasonably 
probable. Further work dealing with these topics could be useful in 
establishing some specific estimates of discharge magnitudes. 

An attempt was made to develop a model of bubble evolution to be 
used in conjunction with ignition and charging theories. Some basic 
relationships giving bubble evolution and movement were developed, 
although time and scope of the study did not allow investigators to 
completely develop a means of predicting bubble size and velocity. 
However, future work in this area would help relate bubble evolution 
to droplet and bubble charging since it appears that breakage of bubbles 
emerging from the surface is a contributing process to static charge 
bui Id-up. 

In conclusion, there are a number of processes that contribute to 
autoignition which appear probable. However, since some of the processes 
cannot be adequately described with the information at hand, much specu- 
lation results. The high energy portion of the critical mass curve comes 
from theory as well as experimentation which makes it quite believable. 
The constants in the equation depend on empirical results and more 
definitive experiments would be desirable. However, experience in many 
tests conducted by Farber indicate that the constants in the equation 
are substantially correct. 
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2.0 INTRODUCTION 



. This study was concerned with the explosion and fire from propel- 

lant tankage ruptures. Such ruptures can be due to equipment failure or 
deliberate action in testing objects to obtain blast yield parameter 

, prediction models for safety studies. The ignition phenomena, wherein 

the studied fuel and oxidizer spontaneously detonate, has been called 
"autoignition". The yields observed in the tests indicated that little 
fuel and oxidizer mixing occurred at the time of ignition. One proposed 
mechanism for ignition at low mixed mass is a spark discharge caused by 
an electrostatic charge which builds up as a result of relative fluid 
motion between fuel and oxidizer. 

2.1 BACKGROUND 

An extensive experimental program, which preceded this study, was 
carried out under the direction of Dr. Farber. Small- and large-scale 
tests were carried out with propellants and many supporting experiments 
were conducted to gain further understanding of the mixing and ignition 
mechanisms; some theoretical development also accompanied the testing. 

The results of the experimental studies are detailed in numerous reports 
and summarized in the final report.* Farber approached the interpreta- 

* tion of the data from two basic directions: physical models and statist!' 

cal models. The physical models led to the Critical Mass Hypothesis (as 
discussed below) while the statistical models were helpful in correlating 
predictions of explosive yields under various conditions. 

The Critical Mass Hypothesis is based on the understanding gained in 
separate experiments on electrostatic build-up and mixing. Stated 
simply, the Critical Mass Hypothesis says that if the two specimens are 
mixed with a specific mixing energy, then a mixed mass will be reached at 
which 100% ignition is probable. In other words, each mixing energy has 
a specific critical mass. The energy released by the boiling of one 
species is associated with a critical mass of about 2300 lb or less. All 
observed explosions should fall on or below the critical mass line. It 
was observed that the S-IVB’ test failure fell just below the critical 
mass line indicating the relationship of Critical Mass Hypothesis to 
failures in actual vehicle configurations. 

The Critical Mass Hypothesis is based on the idea that the relative 
motion of the two fluids, which results from the mixing energy, causes a 

* build-up of a static charge field. This field discharges, causing sparks 
that are sufficient to ignite the mixture. Other ignition sources may be 
present so that the mixture ignites masses that are well below the criti- 
cal mass (as is often observed). No explosion occurs above the critical 
mass since the certainty of autoignition is 100% at the point when the 
critical mass of materials has been mixed. 

t 
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The reality of the presence of static charge build-up was demon- 
strated in Farber's work by the observation of static charge in the 
fluids during mixing. Metal screens placed in containers registered 
measurable potential build-up when propellant constituents were dumped 
together. 

While the concept of autoignition appears to be experimental fact, 
it remains desirable to understand the physical basis for the phenomena 
from the standpoint of first principles. 

2.2 OBJECTIVES OF THiS STUDY 



t . / 


The main objective of this study was to construct a theoretical 
basis for previously observed autoignition phenomena in LOX/LH 2 and 
LOX/RP-1 to whatever extent feasible within the program scope. Theo- 
retical models of competing processes (for ignition and charge build-up) 
were developed to attain the objective. The most probable mechanisms 
for autoignition were identified and described in detail. However, 
development of statistical models for explosion probability and yield 
were not within the scope of this study. 

2 . 3 SUMMARY OF APPROACH 


The study was carried out in four basic areas (with considerable 
overlap between phases). In the first area, the "physical" environ- 
ment* was described as much as possible. The process of bubble formation, 
movement and break-up was modeled as well as the process of solid formation 
and interstage mixing. In the second area, the static charge build-up and 
its relationship to requirements for detonation was described from theoreti' 
cal principles. Other detonation sources (besides electrostatic) were con- 
sidered in the third area and their importance, relative to the electro- 
static mechanism, v/as considered. It is important to stress that so-called 
"external" ignition sources (such as those rising from metal-metal sparks) 
are not considered in autoignition, since such sources are not 
necessarily present. If these sources are present they will cause detona- 
tion at lower mixed volumes than does autoignition, thus giving yields 
below those predicted by autoignition (as has often been observed). In 
the fourth area, an attempt was made to bring all the items from the first 
three areas together to explain the critical mass curve or provide an 
adequate substitute. 


4 


* In this report "physical" environment refers to the various states of 
the fluids that are involved and the dynamics of intermixing. 
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3.0 THE PHYSICAL ENVIRONMENT 


The physical environment governs the mixing and ignition pro- 
cesses. This section describes three aspects of the environment 
that are pertinent to autoignition phenomena: 1) the formation of 

vapor bubbles as the warmer fluid causes the cold fluid to boil, 

2) the physics of mixing the two components with a given energy, 
and 3) the evolution of solid particles as one fluid freezes. 

Bubbles are particularly important since they probably con- 
tribute an essential function to charge distribution and localization 
of discharges. Also, bubbles may be responsible for the vertical charge 
distribution Farber observed with his instruments. 

The Critical Mass Hypothesis is directly related to the dynamics 
of fluid mixing. Fluid mixing and electrostatic charge build-up are 
interrelated in the formulation of the concept of a critical mass. 

Solid particles of oxygen may be the source of oxidizer in the 
liquid hydrogen plug. Hydrogen bubbles associated with oxygen crystals 
could be the detonation mixture in LOX/LH 2 systems. 

3.1 BUBBLE POPULATION DYNAMICS 


Bubble population dynamics involves a description of the bubble 
population in a mass of boiling liquid (e.g., the H 2 plug falling 
into 0?). The formation, movement, and disappearance of bubbles is 
described. 

The details of the model development can be found in Appendix A. 
A general model is described and some pertinent simplified cases are 
further developed. The models are summarized in this section. 

3.1.1 The General Model 


The general model includes five basic processes: 

• Bubble production by nucleate boiling 

• Bubble flow into and out of a region 

• Bubble growth into and out of a size range 

• Bubble break-up 

• Coalescence of two or more bubbles 

The general expression which includes all these factors is: 


G(R) |^+ B(R)n = - V 


• i + s(r,t)f(R) 



n(R')B(R' ^ R)dr 



n(R')n(R")C(R',R" R)dR'dR" 


+ 



- n(R) 



n(R')C(R,R')dR'. 


(3-1) 


where i 
s 

f(R) 

B 

C 

G 


the flow of bubbles (bubble current) 
nucleation probability 
nucleation size distribution 
break-up probability 
coalescence probability 
growth probability 


While this expression includes a general description of all important 
phenomena it does not, in general, provide a convenient means of solving 
for n other than by involved numerical analysis. Specific restrictions 
governing the behavior of the system must be introduced to obtain usable 
models. 


3.1.2 Some Simplifications 

When applying the general equation to specific cases, more concrete 
statements must be made about specific gravity. Frequency functions 
must be detenriined for nucleation, coalescence and break-up; the rela- 
tionship between bubble velocity, size and position must be defined; 
the relationship for bubble growth rate as a function of size and posi- 
tion must be developed; and models for discontinuous re;' ions must be 
matched up. The task is clearly formidable for general models especially 
in obtaining solutions to the equations. However, some cases are rea- 
sonably tractable and provide some insight into the bubble population 
in the mixing propellants. This, in turn can be related to charge 
build-up and discharge across or among bubbles. 

A simplification, common to all the specific models, is the neglect 
of bubble coalescence and break-up. This is reasonable if the bubble 
density is low and no very large bubbles exist. A further simplifica- 
tion occurs if the bubble velocity depends only on size. In this case, 
since there is only one growth function for all regions, it is now 
possible to write the current i as 


i = n V(R). 
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The expression for j(R) can be inserted into Equation 3-1, and 
with the coalescence and breakup terms ignored, the result is 


+ G(R) + V(R) • vn = s(r,t) f(R). (3-2) 

Expressions for G(R), ^(R), S(r,t) and f(R) are needed to make 
the relationship solvable. It then becomes a first order, partial, 
differential equation for n(r,R,t). 

3.1.3 The Two Region Model 

In the case of two liquids coming into contact with one boiling, 
the general equation must be applied to two regions: one where bubble 
evolution and growth occurs in a boundary layer under a supersaturated 
condition, and one where the surrounding liquid is saturated but not 
superheated. 

The Boundary Region 

The growth in the boundary region can be described from empirical 
relationships in the literature (see Appendix A). The radius of the 
bubble as a function of time is 



where K = conductivity 
Py = vapor density 
Q = latent heat 
a = thermal diffusivity 
= superheat 
t = time 

Therefore 

R = K K = constant 

A force balance on the individual bubble provides an equation of 
motion on the assumption that the chief impetus for bubble motion in 
the boundary region is buoyancy: 


I 

1 

4 





4 

t 


A 




A 
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Net Force 


Buoyancy 


Drag 


4 dV 

3 " dt 


7T pg 


6 IT li RV. (3-4) 


The problem becomes one-dimensional in the vertical direction so 
that r becomes Z. Rearranging Equation (3-4) gives 


^ V = S3. 

2p.R^ ^b 


From the growth equation (Equation 3-3) 


dt ■ dR dt ■ dR ^2R J 


(3-5) 


Equation (3-5) can be solved for V(R) so that 


V(R) 




(3-6) 


where 


C = -5-^— 
1 72 

% K 


C = 

^2 72 

K p. 


therefore 


dt 


= V(R) 




(3-7) 


and 


Z = Z. + 




0 \2(C^ + 2) 


(3-8) 


In one dimension. Equation (3-2) becomes 


an 


3t = G(R)|a*V(R)§=S{z.t) f(R) 


an 


(3-9) 
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where 
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6(R) 


dt ■ 2R 


V(R) = 



R^ 


If we assume that nucleation proceeds at a constant rate and a 
single size then 


f(R) = 6(R) dirac delta function 


■» >» 


and 

S(z,t) = S = constant rate of nucleation 
Equation (3-9) can be solved to give 


n(z,R,t) = 


yfTz 


where 


H(t) = 


1 t > 0 

0 t < 0 


r-k(| 


H 

1 

1 

( ji 

\ 

■■ 

/ 

m 


L ' 



(3-10) 


C2K^ 


” ■ 2(C^ + 2) • 

At the interface where Z = h 


n(h,R,t) = — ^ 6 
VAh 


r /k \V4i 



1 

7^1 
>1 = 

1 

H 

‘-(s) 


(3-11) 


For example, if the degree of superheat is 5°C (the real value 
will be greater than zero and less than ZO^C) in liquid hydrogen then 


K - 0.068 m/sec 
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T 


I 


= 3.89 

C 2 = 2.97 X 10^ (tn/sec)"^ 
A = 1 .04 X 10^ m/sec. 


Then a 0.25-in. bubble (such as observed by Farber) will emerge from 
a boundary region of thickness 


h = 



7.9 cm. 


6 3 

Furthermore, if S = 1 x 10 nuclei/sec-m at 2 = 0 (an arbitrary guess) 
then 





i 


_ _ 1 1 .. i«5 0.25 in. bubbles 

n = 1.1 X 10 2 

sec - m"^ 


at the interface. 
The Bulk Fluid 


The following set of equations describes the bulk fluid for the 
simplified case to be considered. Detailed arguments developing these 
equations can be found in Appendix A. 

0 (no mass transfer) 

0 (constant total entropy) 

(vapor is ideal gas) 

where r = gas law constant 

Sy(z) = Cpy In ^ - r In ~ (vapor entropy) 

T^ , P.| reference state 


^ _ 
dt ■ 

^ - 
dt ■ 


Pv(z) = 
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r 




= constant liquid density 
Sl(z) = CpL In 1 ^ 

?2 Reference State 

P(z) = Pg - P|_9Z (hydrostatic pressure) 

These equations are manipulated (in Appendix A) to get the density of 
a bubble 


p(z) = 


Pv(z) Q 


/ 


" TP ,1 

T(z) 1 


^PV P^ ■ ^PL TJ 


where Q = latent heat 

_ ’’o - 
r T(2) 

T(z) = saturation temperature at P(z) 


(3-12) 


T^, P^ = T, P at z = 0 
Since ^ = 0» then M = Mg so that 


R(z) 


.K 

[4ir p(z)J 


With similar arguments for buoyancy forces as in the boundary region 
development 


V 


dl _ 6Tt u R(z) 

dz-- Mg 
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or 


dz V 


(3-13) 


3-7 



where 




and velocity can be obtained from 

V(z) » - f f(z')dz’+ f 5|fi| dz- (3-14) 

•'o "V 

which must be solved numerically since it does not have a closed-form 
solution. After numerical integration R(z), V(z) and, therefore, V(R) 
are known and 


where 


If * F(z) f * V(R) If • 0 


(3-15) 


^ ^ ^ _ y/ V ^ 

dt dz dt ~ ' dz 



Equation (3-15) could then be solved numerically to obtain n(z,R,t) 
in the bubble fluid. The boundary condition for the solution would be 
n(h,R,t) obtained from Equation (3-11). 

3.2 THE SOLID OXYGEN PHASE 


In the hydrogen-oxygen case there is a large temperature difference 
between the two components. If the tanks are at atmospheric pressure, 
the hydrogen will be at 20°K and the oxygen at 90°K. The oxygen will 
tend to boil hydrogen and at the same time will cool rapidly. Since 
the melting point of oxygen is 55°K, which is well above the hydrogen 
temperature, it is likely that solid oxygen will form. 

The heat removed by hydrogen/vaporization is related to oxygen 
cooling and freezing by 

* "c ’0 * "0” "'T) 

where 

. Mj^ = mass of H^ vaporized, 

= mass of 0, frozen, 

0 C R M 

M ' = mass of 0- cooled from T “ to T„ , 

0 d 00 
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Cq = heat capacity of oxygen, 

M " = mass of Oo cooled some other amount aT, 

0 Z 

aT = arbitrary temperature difference (unknown), 
= latent heat of vaporization of hydrogfn, 
x^ = latent heat of fusion of oxygen, 

D 

Tq = boiling point of oxygen, and 
ij* = melting point of oxygen. 


The largest crystal mass of oxygen will be "^ormed #/hen all heat used 
to boil hydrogen is taken from a certain mass of oxygen which is sub- 
sequently frozen (that is = o in the heat balance). 

Therefore, the largest crystal occurs when 


= «o=o (^0® - * «0*0 

, "h ^ -'0 

Chen ^ 




H _ 35 (0.394) + 3.3 
108 


= 0.156. 


The volumetric ratio is 


0.156 (1.41) 
1.2 X 10'^ 


= 183 


and the mole ratio is 


moles hydrogen , 32 ,q , 35 
moles oxygen 2 
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If the hydrogen was in the form of 0.25-in. bubbles, then the 
associated oxygen crystal would have a volume of about 7 x 10"^ cm3 
which, if spherical, would be about 0.09 in. in diameter. Therefore, 
one might envision a 0.09-in. diameter oxygen crystal (or smaller) 
associettd with each 0.25-in. hydrogen bubble. 

Th. oxygen crystals would account for very little of the mixed 
hycrog* I volume, but could be important if not essential to detonation. 
Ir another propellant safety study it was found that if a nitrogen- 
oxyg»n nixture of at least 75% oxygen was mixed with hydrogen, suffi- 
cient solid oxygen was formed to support detonation of the hydrogen 
by a spark. 5 The quantity of 75% oxygen added to 300 g hydrogen was 
75 g wh'ch means that about 16% solid oxygen in the hydrogen is needed. 

Th system could be modeled as follows. The hydrogen mass is a 
plug of fluid entering the oxygen pool. As the plug enters, hydrogen 
b/bbles and solid oxygen crystals are formed and dispersed through the 
hydrogen (because of the very small size of the oxygen crystals). Arcs 
across bubbles or from bubble to bubble must be in the presence of 
oxygen close to the bubbles; the essential oxygen would probably be 
in the form of solid crystals. Thus, it appears that one whould view 
the process as a spark (across bubbles or between bubbles) which vapo- 
rizes some oxygen from the solid and sets off the hydrogen-oxygen gas 
mixt, jre. 

3.3 MODELING OF THE MIXING PROCESS 


When one liquid is released into the other, a mixing process 
begiis. The Critical Mass Hypothesis is based on thp idea that a 
static r.hc’rge is bin'lding up simultaneously with the mixing so that 
a critical mass ’ mixed when sufficient charge is built up to 
guarantee ignition. The critical mass curve is based on the relation- 
ship between mass mixed and mixing energy. The shape of the curve is 
governed by the mass-energy relationship while the specific constants 
ere detetmined from the data on charge build-up. 

The mixing model described here (and detailed in Appendix B) 
shows that the shape of the high energy portion of the critical mass 
curve can be derive from first principles. 

In the hiql energy region it is reasonable to assume that the 
mixing is due tO isotropic turbulence. A relationship between mixing 
energy and quantity mixed can then be developed from Kolmogoroff 
theory o turbulent mixing. 

The ratio of the amount of A mixed with B (= m) and the total 
atrotint of A can be described in terms of time and a constant, which 
uepends on the pow..r density (energy density): 



I 








where 


m _ 
M ■ 




TTJ 

P 



- 1/4 


y = power/unit volume 


Incegration of the relationship over a period of time with a constant 
influx rate dM/dt gives 


-k 


m = k^t 


xr^ (1 - e 


where k^ , k^ are constants. 

Using the relationship for energy 


E = Mgx 

where x = ullage space height and expanding the m(t) relationship gives 

m = KjE^ 

where Kg = constant 

as the first term in the series. It can be shown that other terms in 
the series are small for conditions of interest. This term is identical 
with the high energy term in the critical mass equation. 

Therefore, this theory developed on first principles appears to 
support the empirical interpretation used by Farber to develop the 
high energy term in the critical mass curve. 
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4.0 MECHANISMS OF ELECTROSTATIC IGNITION 


A number of electrical ignition mechanisms can be postulated in 
the rapid mixing of liquid oxygen (LOX) with liquid hydrogen (LH^) or 
liquid hydrocarbon fuels (like RP-1). These include: 

1) electric discharge across combustible fuel-oxygen vapor mixtures 
within the mixture bulk, 

2) liquid dielectric breakdown discharge across both fuel and oxygen 
within the mixture bulk, releasing a combustible vapor mix, and 

3) electric discharge through combustible vapor from or among droplets 
ejected from the surface in boiling. 

A necessary condition for any of these electrical ignition mechanisms 
is a charge separation giving rise to a potential difference sufficient 
to cause a discharge. 

4.1 STREAMING POTENTIALS IN FUEL-LIQUID OXYGEN MIXING 

The streaming potential phenomenon has been long known to exist in 
flows through capillaries and porous plugs, and a quantitative treatment 
has long been available although it does have restrictive assumptions.® 

The phenomenon has been described mathematically and consistency of the 
description has been verified experimentally'^’® for laminar flows of 
fluids through capillaries and porous plugs. 

This streaming potential phenomenon occurs as a consequence of the 
electrical double layer at the interface between two materials. The 
interaction of the two materials results in thin regions which are not 
electrically neutral in the region of the interface- Since the region 
of electrical charge is thin but of finite extent, some of the charge 
gets carried away with the liquid flow at a solid-liquid interface (or 
presumably at a liquid-liquid, liquid-gas, or other interface at which 
there is both close contact and relative motion). The convective flow 
of charge results in the buildup of an electric field which in steady- 
state drives a conductive current of equal magnitude and opposite sense 
to the convective current. This electric field is responsible for the 
streaming potential. 

In the mixing of LOX and either LH 2 or hydrocarbon fuels, we picture 
globules, plugs, or columns of a liquid (Liquid A) plunging into another 
liquid (Liquid B) (and possibly bubbling bac . upward due to buoyancy). 

The convective charge flow at the inl;c rface whore liiere is relative 
motion of the two materials would result in a vertical electric field, 
hence a potential difference between vertically separated points. This 
may be responsible for the potential difference measured by Farber 
between vertically separated screen pairs below the lower liquid surface 
when liquid fuels were poured into LOX and vice versa. Since bubbles 
of the more volatile component (hydrogen in hydrogen-oxygen mixing, oxygen 



in oxygen-hydrocarbon fuel mixing) would be forming most rapidly near 
the interface and would possibly have a combustible admixture of the 
less volatile component, the electric field resulting from the convec- 
tive charge flow could conceivably cause ignition by discharge across 
bubbles. In the following discussion we will try to estimate the 
electric field that could result. 

Consider a plug, column, or filament of Liquid A plunging into 
Liquid B. Suppose Liquid A moves with velocity v and that most of 
Liquid B moves more slowly. The convective current i] through a 
cross section of the column is^ 



where the surface is a cross section of the column and pq is the 
charge density. The charge density satisfies the Poisson equation 
in (Gaussian units): 

7^i|) = -4 ttp /k. (4-2) 

H 

If we take the charge density as resulting from the chemical interaction 
at the interface and maneuver Equations (4-1) ano (4-2), we obtain as 
the dominant term 
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Here, i, is the value of the potential li' at the interface and is known in 
electrochemical and electrokinetic work as the zeta potential. The line 
integral in Equation (4-3) is around the edge of the surface S. 

The traditional expression for streaming potential in a capillary 
involves the Poiseuille expression for laminar viscous flow in a 
cylinder. Although we do not necessarily expect laminar flow in th'* 
dumping of Liquid A into Liquid B, we do expect a boundary layer at the 
interface. We also expect the curl of the velocity at the surface to be 


vxv 
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where scaler v = velocity at the center of the column, 6 = the thickness 
of the boundary layer, and t = a unit vector tangential to the boundary, 
c, of the surface, s. Accordingly, the convective current i] is 


't ' 


-- I 'c 

?ir 
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where 1^ = the length of the perimeter of the flowing column of fluid. 

The convective flow of charge results in a charge separation and 
hence in an electric field that tries to restore charge neutrality by 
conduction. The return conductive flow path could be through both the 
column of Liquid A or through Liquid B. To see what electric field 
magnitudes might result, consider the case where Liquid B is of much 
lower electrical conductivity than Liquid A. If the flow continues 
long enough, a conductive current i 2 will build up so to give a zero 
total current: 




^ 1 , • 0 


Here, conductive current i 2 is given by 



ds - a<E>A, 


where a = conductivity of Liquid A, and 

A = cross sectional area of the flow stream 
<E> = average electrical field across flow stream. 


(4-6) 


(4-7) 
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From Equations (4-5) through (4-7) we deduce an electric field magnitude 


oA64ir 


(4-8) 


This approximate relation is useful for estimating the electric field 
which might develop, although the average electric field could differ 
substantially from the peak value for large area flow columns. For 
large area flow columns, the conductive flow and the electric field 
driving it might be limited to regions near the interface. 
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Dimensional arguments give a boundary layer thickness on the 
order of 

■ ■ fe) ■" 


where u = coefficient of viscosity, 

= mass density, 

V = velocity at points away from the interface, and 

X = distance from leading edge of the surface along which the 
boundary layer develops. 

For our picture of a column of Liquid A penetrating into Liquid B, 
distance x is the distance from the leading end of the column. However, 
Expression 4-9 for the boundary layer thickness should not be taken too 
seriously for x approaching zero nor for 6 larger than the order of one- 
fourth of the transverse flow dimensions, the latter being the laminar 
flow limit. 

We insert Equation (4-9) into Equation (4-8) for a circular 
cylindrical column and substitute plausible values for hydrogen plunging 
into oxygen: 

V = 600 cm/sec (velocity of the plug) 

r, s 0.1 volt = 0.333 X 10 statvoH (estimated zeta potential, 
uncertain by factor -^3 and sign) 

r^ = 3 cm (radius of column) 

X = 5 cm (distance from end of column) 

p = 1.4 X 10'^ gm/cm sec (LH 2 viscosity at boiling point, 20°K) 

K = 1.23 (dielectric constant) 

( 1 / 0 ) = 10^^ ohm cm = 1/9 x 10^ sec (resistivity of liquid hydrogen)* 

2 3 

= 7.1 X 10' gm/cm (liquid hydrogen density). 

The estimated electric field with this set of numbers in Equations (4-8) 
and (4-9) is 

<E> = -3.6 X 10^ statvolt/cm = -1.07 x 10® volt/cm. 


* The estimated hydrogen resistivity is uncertain by orders of magnitude 
due to ionic impurity content. The number used was reported by 
Cassutt^i = loi' ohm cm. 
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Using an electric resistivity 4 orders of magnitude lower would give 
<E> - 1.07 X 10^ volt/cm, which is still a plausible number 

The breakdown fields in hydrogen gas at atmospheric pressure and 
20°K are '^1 x 10^ volts/cm. Thus, breakdown field strengths could 
conceivably be greatly exceeded by the streaming potential mechanism, 
but sharpening of both the model and the physical parameters is needed. 

One interesting property of streaming potentials is that the zeta 
potential, c, can be markedly affected by impurities in the liquid. For 
example, the effect on the zeta potential of adding Th(N03)4 to water 
(shown in Table 4-1) is appreciable for very small quantities of solute. 
In fact, the sign of the zeta potential can actually be changed by 
impurities. This leads to the speculation that reversed polarities in 
Farber's voltage measurements could have been due to build-up of impuri- 
ties in the test system. 


TABLE 4-1 . Effect of Th (N 03)4 ofi Zeta Potential 
Between Water and Glass 


Zeta Potential 
(volts) 

-0.16 

0 

+0.14 


Concentration of Th(N 03)4 
in Water 

(ppm) 

0 

48 

4800 


Marked changes in resistivity and the zeta potential itself could 
greatly alter the charge build-up rates. However, it is not certain 
what kind of impurities one might find in liquid hydrogen or oxygen. 

By the nature of their cold temperatures, the cryogenic liquids resist 
dissolution of most impurities. 

4.2 DISCHARGE ACROSS VAPOR BUBBLES IN THE LIQUID 

Dielectric breakdown within the liquid is plausible for ignition. 
However, it requires not only an electric field which is sufficient to 
cause breakdown, but also requires dissipation of enough energy to vapo- 
rize material in a region of transverse dimensions which are larger than 
the quenching distance (quenching distance is smallest distance from a 
wall or other structure for which a combustion process can continue). 

Walls or liquid materials can "quench" combustion in vapors by robbing 
either heat or reacting ions. The quenching distances for various con- 
centrations of hydrogen in air are given in Table 4-2.^^ A spacing of 
0.06 cm is required for stoichiometric properties with much larger spac- 
ings Approaching 1 cm) for high or low hyoroyen percentages; hydrocarbons 
(like RP-1) would require even greater distances. We anticipate somewhat 
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shorter quenching distances for dielectric surfaces like LOX or LH 2 
than these experimental values measured with metal electrodes, but of 
the same order. See Appendix C for discussion. 


TABLE 4-2 . Parallel Plate Quenching Distances 
for Hydrogen-Air Mixtures 

Hz by Quenching 

Volume, % Distance, cm 


5 

0.60 

10 

0.15 

30 

0.06 

60 

0.20 

70 

0.60 


Liquid dielectric breakdown ignition might result from discharges 
terminating on oxygen crystals in LOX-LH 2 mixing or on wax-state hydro- 
carbon fuels in LOX-liquid hydrocarbon mixing. Further discussion of 
oxygen crystals appears in Section 3.1 under "The Solid Oxygen '-hase". 

The contact potential mechanism gives rise to a finite thickness 
charge layer at the interface between materials, as will be discussed 
further in Section 4.3. Rising bubbles of hydrogen vapor will have a 
tendency to have part of the charge outside them sheared off by the 
relative motion as shown in Figure 4-1. The charge separation results 
in a vertical field component which tends to pull the external charge 
along with the bubbles. This vertical electric field component could 
contribute to the potential difference observed by Farber between 
vertically separated horizontal screens. 

The charge in the liquid just outside a bubble surface is 
approximately 

Ob ' Oj 

where Q« is charge per unit area. Relative motion of bubble and liquid 
attempts to sweep this charge off, resulting in an average convective 
current Jc in the liquid from a dense cloud of moving bubbles given by 

^^c " ”v ''eff % 

where '.’eff is the effective (for carrying charge) flow velocity component 
in the mean flow direction. Here ny is number of bubbles per unit 
volume. I he value of effective velocity Vgff is less than the relative 
velocity of fluid to bubbles, because of boundary layer effects. It is 
a function of both charge layer thickness and flow boundary layer thick- 
ness. An electric field develops whir' tends to keep the charge moving 
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FIGURE 4-1 . Rising Bubble Charge Transport 


with the bubbles. Equating this conductive current magnitude to 
that of the oppositely directed convective current gives 



f 
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where Vav is average relative velocity of fluid to bubbles, s is 
charge layer thickness, and 5 is fluid boundary layer thickness. We 
assume s < 5 . We will take as an approximation for Qa (in cgs 
Gaussian units) 



Where Vq is a potential at the interface that is on the order of the 
difference in work functions for the two materials and ; is relative 
dielectric constant. These approximations give 


u 0 

^av 47 t r6 


We further take the average relative velocity Vgy to be the StOKes 
velocity (valid only for small bubbles). 


V 


av 




9 . 


where m is fluid viscosity, p|_ and pb are liquid and bubble densities, 
and g is the acceleration of gravity. With this Vav we obtain 


vE> 


6 r ^ w 


Inserting arbitrary but plausible parameters, 

0 = 0.009 sec'^ = lO"^^ (ohm cm) ^ 
a = 0.1 
f. = 1.46 
( r /:0 = 10 . 

3 

= 1.14 gm/cm 

Po = negligible 

“ 2 
g = 980 cm/ sec 

Vq = ^ statvolt = 0.1 volt 

-3 -1 

'I = 9 X 10 gm (cm-sec) 
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we obtain 


= 3.6 X 10^ statvolts/cm = 1.1 x 10^ volts/cm 

This field strength compares witii a dielectric breakdown field of 
10^ volt/cm in 20'^K hydrogen gas. This suggests that discharge 
across rising hydrogen bubbles is possible by this mechanism. 

4.3 ELECTROSTATIC CHARGING EFFECTS AT SURFACES OF MIXING LIQUIDS 

If two liquids whose vapors from combustible mixtures are vio- 
lently mixed (as by dumping one into the other), the cover gas-liquid 
interface seems a likely place for ignition. For volatile materials, 
a combustible '■ixture at many locations over the surface is probable. 

One is therefc.e led to seek a mechanism for electrical discharge of 
sufficient intensity for ignition at the surface. 

Experimental work of Harper'"" verified charging of liquid droplets 
formed at a number of liquid surfaces. His data was obtained by bub- 
bling nitrogen gas up through the liquids. Charged droplets of both 
signs were seen (though not in equal numbers) when froth bubbles broke. 
While the exact mechanism for this charging cannot be deduced from 
these experiments, the experiments do seem to prove that droplets 
charging so formed does occur. We believe that contact potential 
mechanism is responsible for the charging observed by Harper. If an 
electric field is present, an induction process can also contribute to 
droplet charging. Harper, however, attributed his droplet charging to 
shear in electrical do’ble layers. 

qualitatively, a space charge layer is known to occur near the inter- 
face between dissimilar materials. The space charge region supplies the 
electrical potential difference needed to give zero net current at equi- 
librium between dissimilar materials in contact. Treating the case of 
metal -semi conductor contacts, Schottky proposed a model in which the 
width of this space charge region is proportional to the square root of 
the contact potential and inversely proportional to the square root of 
the number density of donor (acceptor) atoms for an n (p) type semiconduc- 
tor. If this qualitative dependence on donor atom density holds, then the 
space charge region becomes thick for highly resistive, very pure cryogenic 
liquids. A thick region of space charge in the liquid oxygen around a 
hydrogen bubble when it breaks the surface gives an opportunity for 
charged oxygen droplet formation. Presumably there would be opposite 
sign charging on high mobility hydrogen molecules. 

This postulated liquid oxygen droplet charging at the surface of 
the LOX-hydrogen mixture could conceivably lead to ignition in a number 
of ways: 
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1) Electrical discharge between an oxygen droplet and the bulk 
liquid as the droplet is ejected away in froth bubble breakage, 

2) Electrical discharge between a charged droplet and the bulk 
liquid as the charged droplet falls back into the liquid, 

3) Discharge between flying droplets with dissimilar levels of 
charging, 

4) Miniature "lightning" flashes if cover gas flow over the surface 
or difference in mobilities separates charged oxygen droplets 
from oppositely charged hydrogen gas, 

5) Droplet-droplet discharge by the following sequence of events: 
Cover gas flow over the LOX-hydrogen mix, difference in 
mobilities, or gravitational effects separate oxygen droplets 
from oppositely charged hydrogen gas. Droplets formed under 
the non-space-charge-neutralized oxygen droplet cloud get 
opposite sign charge by induction. Droplets of opposite sign 
spark discharge in passing. 

A sketch showing droplet charging at the liquid surface is shown 

in Figure 4-2. 


HYDROGEN BUBBLE 
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FIGURE 4-2. Surface Froth Droplet Charging Effects 
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For mixing LOX and RP-1 , similar postulates may still explain 
surface autoignition. Oxygen vapor bubbles rise rapidly to the 
surface with contact potential charge distributions in the RP-1 
layers around it and in the oxygen vapor of the bubble. The RP-1 
droplets get a charge when the froth breaks at the surface. Some 
of the RP-1 around a rising oxygen bubble may be wax, just as some 
of the oxygen around a risirtg hydrogen bubble may be crystalline. 

It should be noted that in the autoignition work of E. A. Farber 
et al., the ignition source location in Figure 80 of Reference 1, 
suggests a surface ignition phenomenon. 

4.3.1 Contact Potential Effects in Cryogenic Fluid Mixing 

The conduction mechanism in liquid hydrogen and liquid oxygen 
probably involves impurities, but it is not clear (to this writer) 
whether the charge carriers are ions or electrons or holes. We favor 
electrons or holes as the contributors, since ion mobility should be 
low at cryogenic temperatures. Also, at cryogenic temperatures it 
seems likely that there is a fair amount of spatial order of molecules, 
i.e., on a scale of several to several hundred intermolecular distances. 
Hence, it seems appropriate to use the terminology and models of solid 
state semiconductor physics to these liquids (e.g., conduction and 
valence bands, conduction band electrons and valence band holes, donor 
and acceptor impurities, etc). At the same time, it seems appropriate 
to regard conduction band electron wave functions (if this is a valid 
concept) as localized over a distance on the order of the molecular 
order range. Thus, it is appropriate to speak of a spatial variation 
of electron number density on distances of this order; we emphasize 
that these conjectures are (at least to us) speculative. 

The high resistivity of cryogenic liquid hydrogen'^ of high 
purity suggests that the number density of charge carriers is very 
low. If conduction is by electrons or holes, the resistivity should 
be extremely sensitive to donor or acceptor impurities, because the 
number density of intrinsic charge carriers is so low. Measurements 
of L. Cassutt, et al.,'^ showed insensitivity of resistivity to air or 
nitrogen impurities, which suggests only that the air or nitrogen 
impurities were not good donors or acceptors in liquid hydrogen, i.e., 
they did not ionize. We did not uncover corresponding information for 
liquid oxygen. 

The charge density near the interface needed to produce the electro- 
static potential difference obeys Poisson's equation (written here in 
cgs Gaussian units) : 
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Here, the assumed spatial dependence of the number density n^ of 
charges of species i is given by: 


"(W ' "io 


q^V(r) 


(4-12) 


where V is electrostatic potential, q^- is charge per particle of species 
i, T is absolute temperature, and e is the dielectric constant of the 
medium. This relationship can be considered a classical Boltzmann 
statistical dependence for ions or as a low concentration limit for 
electrons and donors or holes and acceptors, as we shall discuss later. 

To see the nature of the charge distribution near the interface, it 
will be sufficient to consider the region x > 0 when the interface is 
at X = 0 and to assume that for x large, V, p, and dV/dx approach zero. 
With only an x dependence of electrostatic potential, V, and these 
assumptions. Equation (4-1 1) yields a first integral 

If we take only two charge species of positive charge q and negative 
charge -q and take the positive x side as medium 2 when V] - V 2 ^ o (so 
dV/dx < 0 ), then Equation (4-13) gives 


dx 


^8TTkTnQy/2 




(4-14) 


where no is the number density of either sign charge at large x. 
Separation of variables in Equation (4-14) and integration yields 


[exp (qV/kT) + l][exp (qV /kT) - 1] 

log 

[exp (qV/kT) - l][exp (qV^/kT) + 1] 


<x 


where 


^ ckT 


Here, V is the value of V at the plane of the interface. 
0 


(4-15) 


(4-16) 
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Similar treatment of the other side of the interface gives 




[exp (qV/kT) - l][exp (qV /kT) - 1] 

log ; ^ - KX (4-17) 

[exp (qV/kT) - l][exp (qV^/kT) + 1] 


where 

V = V - (V^ - V^) 

Vq - V (x = 0) 

/S.q^n \ 1/2 (4-18) 

\ekT y 

n^ = Number density of current carriers for X '^< o 
: = Dielectric constant for x < o. 


Note that V] - V 2 will be determined at equilibrium by Equation (D-5) 
in Appendix D, "Contact Potential Fundamentals" and will be essentially 
a work function difference. Vq will be between zero and - V 2 , and 
can be calculated by requiring continuity of dV/dx across the interface. 

The space charge region will be of thickness 


{<) 


-1 


ekT 
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in Medium 1 (x > 0 ) and of thickness 


= (<)'^ 


(4-20) 


in Medium 1. To understand this, consider the case V -- kT. Then 
Equation (4-15) gives 


V(x) e 


(4-21) 
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A graph of qV/kT versus <x according to Equation (4-17) for several 
values of qVo/kT can be found in H. R. Kruyt, Reference 9. 

Note that Equation (4-19) implies thick space charge boundary 
layers for interfaces at a low temperature in a media with a low number 
density of mobile current carriers. It is necessary to show that the 
number density of charge carriers can be sufficiently low that the 
space charge region can be on the order of the thickness of the liquid 
oxygen film when a hydrogen vapor bubble emerges at the surface. We 
would then expect significant charging of droplets formed when the film 
breaks. If the space charge region is only a few molecule separations 
thick, it would appear that charge separation in the film breakage 
would be unlikely. 

We can calculate the charge per unit area in Medium 2 at the 
interface from 




(4-22) 
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For |qV^/kT| 1, this gives 
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The charge per unit area is larger when the spatial separation is small. 
This was to be expected, since it is the dipole moment per unit area 
which gives the electrical potential step to satisfy the condition for 
no net electron flow between media [Equation (D-5), Appendix D]. 

Suppose that a film of thickness Sf of liquid oxygen ovt ' a hydro- 
gen bubble ruptures. Suppose further that some area A of this film then 
forms into a droplet of radius r, carrying along the change q = iA from 
the piece of film. Assuming volume conservation, we have 
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The radial electric field at the droplet surface should have a magnitude 
given by 




(4-27) 


Us ng the Approximation (4-25) for ; gives 


E = 



“ ' IV h ' 


(4-28) 


Suppose a film area A = 10^ Sf^ of oxygen forms a droplet. If 
= 1.507 for oxygen, the space charge layer thickness Sp is 10"5 cm and 
Vq is 0.05 volt, then the electric field at the droplet surface from 
Equation (4-28) would be 22,500 volts/cm. Parallel plate breakdown 
voltage in air at standard conditions is about 31,500 volts/cm and in 
hydrogen gas is 17,500 volts/cm; this droplet charging to the verge of 
electric discharge is quite plausible. 

We note that at atmospheric pressure oxygen melting temperature, 

T = 55"K, Equation (4-19) gives a charge-layer thickness 


$2 = 48.7 cm (4-29) 

where current car'ier density no is in cm"3. For a few charge carrier 
densities this gives 


4.87 

x 

10-® 

cm 

for 

irt+14 ^„-3 
n^ = 10 cm 

0 

4.87 

X 

1 

o 

cm 

for 

n^ = lO'^^^ cm'^ 

4.87 

X 

10 

cm 

for 

n = lO"^® cm’^ 
0 


(4-30) 




4-15 


For comparison with semiconductors, room temperature (300'’K) germanium 
and silicon have intrinsic (as opposed to impurity) electron carrier 
density no = 2.5 x 10l3 cm-J and 1.4 x 10^0 cm-3, respectively.^^ 

Following our postulate of doped semiconductor behavior, let us 
see what impurity carrier density may occur. Impurities of interest 
may be of the donor, acceptor, trap type.^^ Donor impurities have 
electrons at energies slightly below the conduction band, and readily 
supply electrons to the conduction band by ionization to form a posi- 
tive ion. Acceptor impurities can capture electrons from the valence 
band to form negative ions, leaving vacancies or "holes" in the 
valence band as current carriers. Traps can capture or lose electrons 
also, altering the carrier density. 

Kittel^‘^ shows that a perfectly pure nonconductor has an intrinsic 
carrier density ngv’’) given by 

n^('^ = 2 (2TtkT/h^)^/^ [-Eg/(2kT)] (4-31) 

where mg and m^ are effective masses of conduction electrons and valence 
band holes, respectively. Eg is the energy gap width, T is absolute 
temperature, and h is Planck’s constant. If the material is doped with 
donor type impurities only, Kittel shows that the electron carrier density 
will be approximately 

= [2 (2TrmgkT/h^)^''^ exp [-E^/(2kT)] (4-32) 

where is number density of donor impurities and E^ is the minimal 
energy increment to move the donor's electron to the conduction band. 

Since we do not have good information on electron energy levels 
and band structure in liquid oxygen and liquid hydrogen, we shall only 
say that plausible values of E_, E^j, N^, mg, and m^ can give values of 
Oq in the range used in Examples (4-30). 

A charge distribution in a leaky dielectric of relative dielectric 
constant e and conductivity o relaxes with a time constant t given (in 
MKS units) by 


T = (4-33) 

a 

IP p .1 

where Eq is >'8.85 x 10 coulomb joule meter . If the value of o 
for oxygen is on the order of lO'l' (ohm cm)"' or 10"’^ (ohm meter)"' 
as reported by Cassutt et al . , for pure liquid hydrogen, then the 
relaxation time in oxygen would be 1.3 x 10® sec. This says that the 
charge separation would probably proceed only a very short way toward 
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the equilibrium distribution at the interface, but very little 
separated charge could coalesce in film rupture. If the impurity 
level was such as to make o = 10“8 (ohm cm)-l, then the time constant 
would be 0.13 seconds, less than typical bubble rise time. For ref- 
erence, the conductivity of good conductors is ''I 06 (ohm cm)"l, that 
of good insulators is 10"^^ to 10~22 (ohm-cm)-l. Solid state semicon- 
ductors have conductivities of order lO"^ to 10 ^ (ohm cm)‘l. 

The characteristics which go together are 


high electrical conductivity, 

thin regions of charge at 
interfaces, 

^ high charge density at 
interfaces, 

fast charge relaxation, 


low electrical conductivity, 

thick regions of charge at 
interfaces, 

low charge density at 
interfaces, 

slow charge relaxation. 


It appears that the parameter domain for hydrogen-oxygen froths in large 
scale violent mixing may be such as to give droplet charging, but we 
have many gaps in our information. Also, the interface charge distribu- 
tion does not need to have proceeded very far toward equilibrium to give 
some droplet charging. 


4.3.2 Droplet Discharge Possi b ilities 

In Equation (4-28) and the discussion following it, supported by 
Examples 4-30, it was shown as plausible that a droplet be charged to 
near electrical breakdown in film breakage. A much lower level of 
single droplet charging could produce breakdown from the collective 
action of all their charges, similar to the natural phenomenon of 
1 ightning. 


Consider a region of thickness h over the liquid surface carrying 
ny droplets per unit volume with charge q^ on each droplet. Then the 
electric field at the edge of this layer from droplets inside is (in 
rationalized MKS) 


Qjh 

v^ d 
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Suppose that each droplet is charged so that its surface field is some 
fraction ■ of breakdown field En: 
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Further assume number of droplets per unit volume is such that 
droplets occupy fraction B of the total volume in the layer: 

"v I ^*^d^ " ® 

If we insert Assumptions (4-35) and (4-36) into Equation (4-34), we 
can solve for the thickness h needed to give electrical breakdown 
field at the layer edge, namely 


h = . (4-37) 

Hence, if droplet surface field is fraction c = 0.1 of breakdown and 
the droplets occupy fraction g = 0.01 of the space, and droplets a>"e 
10"2 cm in radius, a layer h = 6.7 cm thick of these suspended droplets 
would give electric breakdown field at the surface. 

It is obvious that an electric field strength appro'-'-'-in' 
magnitude would contribute to inductive charging of ne.v o . 's 
with opposite sign as they are formed at the liquid surface. While 
ejection of droplets of opposite charge would tend to reduce the 
electric field, it also would give a good opportunity for droplet- 
droplet discharge. 
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5.0 ALTERNATIVE IGNITION SOURCES 


5.1 COSMIC RAY IGNITION MECHANISMS 


High energy charged particles leave a path of ionized atoms and 
recoiling electrons and atoms along their trajectories. The heating 
and ionization would be conducive to ignition if the energy deposit 
per unit particle track length were sufficient. Since the energy 
deposit per unit track length is orders of magnitude smaller than that 
required for ignition from a discharge, the cosmic ray ignition can be 
ruled out as highly improbable. 

In an experiment on hydrogen ignition, Yates et al.,^® varied 
spacing between small electrode pairs to determine minimum ignition 
energy and quenching distance. They determined that the energy needed 
in a discharge for ignition decreases as the electrodes are brought 
closer together until the quenching distance is reached. Further, 
they found that at electrode spacings that are less than the quenching 
distance, the minimum ignition energy goes up because of energy and ion 
losses to electrodes. With this in mind, one suspects that if the 
energy deposit per unit track length from a charged particle in a com- 
bustible mixture exceeds (smallest discharge energy)/(electrode spacing) 
needed for ignition, the mixture will ignite. The condition for ignition 
is then 

ds ^ 

where is minimum energy for ignition with electrode separation l, or 

dM 

where dE/dM = energy deposit per unit mass per unit area traversed and 
p = the mass density of the medium traversed. 


The minimum ignition energies and corresponding electrode spacings 
from Yates et al. for a 65% hydrogen, 35% oxygen mixture are shov/n in 
Table 5-1. 
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T ABLE 5-1 . Minimum Ignition Energy and Corresponding Electrode 
Spacing for 16°C 65% Hydrogen, 35% Oxygen Mixing at 
14.75 psia 


Electrode Ignition 

Spacing Energy 



(cm) (joules) 


Ignition 

Energy/Length 


(joules/cm) 


Ignition Energy/ 
(Mass/area) 

iig/Mi 

[joules/(gm/cm^)2 ] 


1.016 

3.5 

X 

10-3 

3.4 

X 

10’^ 

6.5 

0.508 

0.75 

X 

10'3 

1.5 

X 

10"^ 

2.8 

0.330 

0.7 

X 

O 

1 

OJ 

2.1 

X 

10-3 

4.0 


We can compare the values of Ejg/(p)i) with the maximum values of 
dE/dM along tracks of selected charged particles in the hydrogen-oxygen 
mixture: 


dM 


* ’O'’'’ 


for *H ions at 0.25 MeV 



7.764 X 10’^ 

gm/cm^ 


for ions at 39 MeV. 


Thus, we see that for charged particle traverses we have 


dE 

dM 


<< 



and the mixture will not ignite. Therefore, cosmic ray ignition appears 
improbable, unless the experimental minimal ignition energies for electric 
discharges reflect son:s controlling condition not present for cosmic rays. 

5.2 CHEMICAL REAC'^ONS 


One major possible source of heat would be the rapid conversion of 
ortho to para hydrogen which liberates about 339 calories per mole 
(cal/mole) .’® At cryogenic temperatures the conversion is catalyzed by 
the presence of a naramagnetic substance such as oxygen.*’ This liberated 
heat has long been known as a nuisance in liquid hydrogen st.i,age since 
it causes rapid evaporation. Considering this possibility, a check was 
made which revealed that propellant hydrogen is not-mally catalyzed to 95% 
para to alleviate such evaporation problems. Since the equilibrium composi- 
tion would be about 99.9% para at 21“K, about 4% of ♦'he hydrogen would 


convert on contact with the oxygen providing about 13.6 cal/mole (or 
57 joules/mole) of hydrogen added. About 13 joules/mole is needed so 
the mechanism is plausible, but not probable, since large heat losses 
from evaporation will prevent heat build-up to ignition temperatures. 
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6.0 CONCLUSIONS AND RECOMMENDATIONS 


Partial explanation and substantiation of autoignition (including 
the Critical Mass Hypothesis) has been possible. Some contributing 
and alternative processes have been identified as mechanisms of auto- 
ignition; electrostatic ignition stands out as the principal mechanism. 

It is likely that a number of electrostatic processes are involved. 
Although probability values for various processes were not calculated, 
order of magnitude estimates were made which indicated reasonable proba- 
bility on a qualitative basis. Some questions remain unanswered and 
some additional work is recommended if a better understanding of the 
process is needed. 

6.1 EXPERIhCNTAL RESULTS EXPLAINED THEORY IN THIS STUDY 

The high energy term of the critical mass equation was determined 
empirically^ to have thp form: 

m = K 2 E^, 

where m = mass mixed, E = dimensionless mixing energy (E = mixing energy/ 
boiling energy) and K 2 = a constant. A mixing theory developed in this 
study indicated the same form of expression. The "ow energy term. 


m = ^ , 

a semi -empirical term based in part on the oscillating plug concept,^ 
was not derived in this study. Since the low energy term is only 
important below the boiling energy, it is not considered very important 
to real autoignition situations because the mixing energy will normally 
be much higher than the boiling energy. No tests have been conducted 
below the boiling energy (see Figure 13, Reference 1). 

Predictions of streaming potential charging, based on the bobbing 
plug concept, indicate that sufficient field can be induced to cause 
ignition. Furthermore, the charging of droplets at the surface accom- 
panied by droplet to dr-oplet discharge is probably a contributing 
mechanism (although lack of specific data make this somewhat specula- 
tive). Some substantiation of this is found in observations of ignition 
points on or near the surface in some tests. 

6 . 2 EXPERIMENTAL RESULTS NOT EXPLAINED IN THIS STUDY 

Cne area that is greatly lacking in data is fundamental measurements 
of charge build-up as a function of mixing. The value of constants in 
the critical mass equation is determined largely from the LN 2 /RP-I mixing 
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tests which indicated a critical mass of 2300 lb at boiling energy for 
ignition of LOX/LH 2 and 2800 lb at boiling energy for LOX/RP-1 . The 
similarity between this data and large scale "real" situations is not 
readily apparent and cannot be readily supported by theory (nor dis- 
proved). Empirical support, however, is reasonably good considering 
that no test has yielded explosions above the established critical 
mass curve. Thus, the results are plausible but theoretical explana- 
tion is not in hand. 

Assessment of probabilities associated with various mechanisms was 
outside the scope of this study. While some mechanisms appear to have 
reasonably large probability, assessment of others has not been 
feasible. Most of the mechanisms which showed reasonable probability 
in LOX/LH 2 are improbable for LOX/RP-1 mainly because of the much larger 
ignition requirements. Nevertheless, the S-IVB^ and PYRO 25,000 lb 
LOX/RP-1 experiments show definite evidence of the reality of bome 
mechanisms. Droplet-droplet charging may be one possible mechanism 
(see discussion on droplet-droplet discharges in Section 4.3) but 
this is very speculative. 

6.3 RECOMMENDATION FOR FURTHER WORK 


If a more complete theoretical understanding is a highly desired 
goal we would recommend additional theoretical and experimental work. 

Whether the explanations are now sufficient is a decision best left 
to those who wish to apply autoignition principles to real applications. 

If more complete explanation is needed then the work described below 
would be highly useful. 

6.3.1 Experimental Work 

More data on mixing energy versus charge build-up is needed. The 
LN 2 /RP-I experiments should be expanded and additional tests with two 
cryogenics (perhaps LN 2 /LH 2 or LN 2 /LOX) should be carried out. The 
position of the critical mass curve (as opposed to its shape) is 
inherently dependent on the results of such tests. It is conceivable 
that the initial mass could be much higher and that external spark 
sources are always sufficient Lo cause premature ignition. 

Work is also needed to characterize impurities in the propellants 
and their role in autoignition. In many cases we have noted that electro- 
static properties (e.g., zeta potential, resistivity) are very sensitive 
to the level of impurities in the substances considered. Autoignition 
could conceivably be greatly enhanced or inhibited by the presence of 
impurities. 

6.3.2 Theoretical Work 

Further experimental work could be accompanied by theoretical work 
to interpret data. In addition to this, expansion of theoretical work 
carried out in this study would be very useful, in some cases. Specifically, 
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further work with bubble evaluation models and with prediction of 
charging effects is recommended. Expansion of these efforts was 
started in this study but time and funds were very limiting. It is 
important to note that expansion of these areas would not be very 
productive unless carried out on a fairly large scale (1 to 2 man 
years of effort), since some break-throughs, not easily acquired, 
would be needed. 


APPENDIX 


BUBBLE POPULATION MODEL 


A.l GENERAL MODEL 


A complete description of the bubble population in a liquid is provided 
by specifying the distribution of bubbles according to their position and 
size at any time. In order to do this mathematically, let d^r denote a 
volume which is small compared to the dimensions of the liquid container, 
but large enough to contain many bubbles. Then, define a function n(r,R,t) 
such that 


I average number of bubbles with 
n(r,R,t)d^r dR = < center located in d^r, and radius 
(lying in dR about R, at time t. 

In order to derive an equation for n(r,R,t), suppose that n(r,R,t) is known 
for some time t, and compute its value at a later time t + At. Clearly, 

n(r,R,t + At)d^r - n(r,R,t)d^r dR = net change in n (A-1) 

during At. 

The net change in n is due to several processes, including 

1) production of bubbles in d^r dR by nucleate boiling, 

2) net flow of bubbles having a radius R in the range dR into d^r, either 
by motion of bubbles relative tc the liquid or by motion of the liquid 
itself, 

3) growth of smaller bubbles into dR and growth of bubbles in dR into 
larger bubbles, 

4) breakup of bubbles in dR into two (or more) smaller components, and 

5) coalescence of two (or more) smaller bubbles to form one of radius R 
in the range dR, and coalescence of bubbles of radius R with one (or 
more) others (s) to form a larger bubble. 

An exact description of all of these processes is probably not within reach 
of existing theory. However, formal expressions may be written for all of 
them, as will be done here: 

1) Let s(r,t)d^r denote the rate that bubbles are produced in d^r at 
time t. Let f(R)dR denote the probability that a given bubble is 
produced by nucleation with radius R in the range dR. Since f(R) is 
a probability distribution, it is normalized so that 
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/ f(R)dR = 1. 

0 

Ip. terms of s and f, the number of bubbles produced in d^r dR during 
a (small) time At is s(r ,t)f (R)d3r dR At. 

Let j^(r,R,t) denote the current of bubbles. Then, by familiar argu- 
ments, the net flow of bubbles of radius R in the range dR into d^r 
during time At is given by - V*j_(r,R,t)d3r dR At. 


Let G(R -* R',At)dR' denote the probability that a bubble of radius R 
will grow to a radius R' in the range dR' during At. Then the number 
of bubbles in d^r with radius R in the range dR, which grow to a new 
radius R' in the range dR' during At, is given by n(r,R,t)d3r dR • 

G(R ->• R',At)dR'. Thus, the total number of bubbles leaving dR by 
this mechanism during time At is given by 


n(r,R,t)d^r dR 


•/ G(R 
R 


R', At)dR'. 


By a similar argument, the total number of bubbles entering dR, due 
to growth from R' < R, is given by 

J n(r,R' ,t)G(R' ^ R, At)dR' • d^r dR. 

0 


To proceed further with these expressions, a model is needed to 
describe G(R' R,At). The simplest model is a deterministic growth 
model which is based on the assumption that the growth rate of the 
bubble radius depends only on its present radius, i.e., 

^ = G(R) 

for some function G(R). For this model, R(t + At) must have the 
value R(t) + G[R(t)] • At, i.e., G(R ^ R',At) = 6[R' - R - G(R)At], 
With this expression for 6(R -> R',At), the above integrals are easily 
evaluated to give 

growth out of dR = n(r,R,t)d3r dR, and 
growth into dR = n(r,R - G(R)At,t)d^r dR. 


Therefore, since At is small 

growth into dR |^n(r,R,t) - G(R) ^ (r,R,T)AtJ d^r dR. 


In the expression for bubble breakup, we include only the spontaneous 
breakup rate which is a function of R alone. (At higher bubble 
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densities, stimulated breakup, which is caused by the passage of one 
nr more bubbles near another, may also contribute, but this process 
will not be considered here.) We let B(R ->• R',R")dR'dr" denote the 
rate at which bubbles of radius R break up to form a pair of bubbles 
with radii R' in the range dR' and R" in the range dR". Then, the 
number of bubbles leaving d^r dR during At, due to breakup, is given 
by n(r,R,t)d3r dR • B(R) • At, where 

R R 

B(R) = // B(R ^ R',R")dR'dR". 

0 0 

dR during At, due to breakup 


■ d^r • At 


B(R* R".R)dR''. 


5) Coalescence is basically the inverse of breakup. Here, two bubbles 

of radii R' and R" merge to form a new bubble of radius R. (At higher ; 

bubble densities, interactions between three and more bubbles may ■ 

also contribute, but this process will not be considered here.) Let 

C(R',R" R)dR denote the rate at which bubbles of radii R' and R", 

both located in d^r, coalesce to form a new bubble of radius R in 

the range dR. Then, the number of bubbles leaving dR during At, 

due to coalescence, is given by } 


Similarly, the number of bubbles entering 
of larger bubbles, is given by 

/ n(r,R',t)B(R' R)dR' 
R 


where B(R' R) = /“b(R‘ ^ R,R")dR*' + T 
R 


n(r,R,t)d^r 



n(r,R' ,t)C(R,R' )dR' 


where 


At 



C(R,R' ^ R")dR". 


Similarly, the number of bubbles added to dR, due to coalescence of 
two smaller bubbles during At, is 

n(r,R' ,t)n(r,R",t)C(R',R" ^ R)dR'dR" • dR • At. 


At this point, formal expressions have been given for each process 
that contributes to the right hand side of Equation (A-1). When these 
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expressions are substituted, each term in he resulting equation contains 
the factors d^r dR which can be cancelled. Then, if the equation is 
divided by At and the limit At -*■ o is taken, the result is 

It ^ ^ ‘ i + s(r,t)f(R) + n(R')B(R' R)dr' 

R 


n(R')n(R")C(R',R" ->■ R)dR'dR" 


- n(R) 


This equation is not yet tractable, since it is necessary to specify the 
current i separately, and to provide models for G(R), B(R' T}, 

C(R',R" -► R) and f(R). Thus, some simplifications to this general 
equation, which lead to more tractable forms, will be discussed in the 
following subsection. 

A. 2 S IMPLIFICATIONS TO THE GENERAL MODEL 

At low bubble densities the coalescence terms will be negligible 
compared to other terms, which lead to 

91 ^ = - V • i +.s(r,t)f(R) + n(R' )B(R' R)dr’. 

■^R 

If, in addition, no very large bubbles are present in the system, the 
spontaneous breakup terms will be small, which lead to 

It If " V • i + s(r,t)f(R). 


n(R')C(R,R*)dR‘. 


(A-2) 



If the bubble velocity depends only on size, the current i can be 
written as 

j(r,R,t) = v(R)n(r,R,t). 


In this case, the equation for n reduces to 

It ^ If T(R) ■ = s(r,t)f(R). (A-3) 
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in this form the equation is tractable if models are available to specify 
G(R), v(R) and f(R). In fact, since 

G(R) = ^ and v(R) = ^ 

for a bubble of radius R and position r. Equation (A-3) simply states that 

s(r,t)f(R). (A-4) 

Here, ^ denotes the total time derivative along the characteristic curves, 
R = R(t) and r = r(t), which are determined by integrating 


^ = G[R(t)] 

(A-5a) 

v[R(t)]. 

(A-5b) 


The solution of Equations (A-4), (A-5a) and (A-5b) for a particular case 
is given in Section A. 3. 

A. 3 SOLUTION FOR TWO ADJACENT FLUIDS - GENERAL 


Suppose that two fluids of different temperature are brought together 
at time t = 0. Idealize the geometry so that Fluid 1 occupies part of the 
upper half-space 0 < Z < h, while Fluid 2 occupies part of the lower half- 
space -b < Z < 0. Suppose the temperature of Fluid 1 is less than that of 
Fluid 2 and that both fluids are initially saturated liquids. A thermal 
boundary layer will form in both fluids with nucleate boiling occurring in 
the boundary layer of Fluid 1 and (possibly) ice crystals forming in the 
boundary layer of Fluid 2. Due to buoyancy, vapor bubbles of Fluid 1 will 
rise out of the boundary layer as saturated vapor and into the bulk of 
Fluid 1. IT Solid 2 has greater density than Liquid 2, any ice crystals 
formed in the boundary layer of Fluid 2 will sink into the bulk of Fluid 2 
as saturated ice (Figure A-1). 

The description of the bubble population requires consideration of 
two regimes: 1) bubble nucleation and growth in the boundary layer, and 
2) bubble rising and growth above the boundary layer. These two regimes 
require separate analyses which are outlined below. 
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FIGURE A-1. A Two-Fluid Model 


A. 3.1 The Boundary Layer 

In the thermal boundary layer, which forms when the two fluids at 
their boiling points come into contact, the colder fluid will be super- 
heated. The formation and growth of bubbles in a superheated region have 
been studied in some detail” with the conclusion that the radius of a 
growing bubble expands as 

R(t) = Kt^^^ (A-6) 

where K is a constant depending on fluid properties and the degree of 
superheating. With the expression for R(t) known explicitly, some of 
the factors appearing on the equation for n (r,R,t) can be determined. 

In particular, from Equation (A-6) 

^ - K .-1/2 _ 
dt ■ 2 ^ = 2R ’ 

Thus, from Equation (A-5a), we have 
G(R) . If • 

The velocity v(R) can also be determined if we assume that the bubble 
velocity is vertical and is due to the combined effect of buoyancy and 
gravity forces. Then, 
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(A-8) 
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(A-9) 


net 

force 


buoyancy drag 


where 


Thus, 


P|^ = bubble density 
p = fluid density, and 
y = viscosity. 


^ + -ilL_ V = £3. 

2pjjR^ Pb 


(A-10) 


Dividing Equation (A-10) by Equation (A-7) gives 

^ ^ V = C R 

dR R '' 


where 


C, = -p- and Cy = 

K Pfa K pjj 


(A-11) 


Here, we take p^ = constant throughout the boundary layer since the boundary 
layer is thin. The solution of Equation (A-11) is then 


v(R) = 




(A-12) 


which provides v(R) (assumed in +Z direction) for Equation (A-4), 
Equations (A-5bT, (A-12) and (A-6) we also obtain 


From 


dt 


C^K‘- 


+ 2J 


(A-13) 


where Z is the vertical component of r. Thus, 
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M 


(A-14) 


C IT 

= ^0 + ?(C“+ T[ 

where Zq is the height at which the bubble first forms by nucleation. 

As an illustration of how Equation (A-4) can be solved for n(Z,R,t), 
we consider the case where one bubble nucleates at Z = 0 at time t = 0. 
Then, according to Equation (A-4) 

^ = 0 for Z > 0, t > 0. 


Thus, n is constant along the curves given by solving Equations (A-5a) 
and (A-5b). Since these have already been solved [Equations (A-6) and 
(A-14) above] for the present simple model, we have for one bubble 


n(Z,R,t) 


2(C^ + 2 ) 


^ 6^R - 


(A-15) 


where 6 = a dirac delta function. 


If, instead of one bubble, we consider the continuous production of bubbles 
at Z = 0, starting at time t = 0 at a rate of S bubbles per cm3 sec, then 
we integrate Equation (A-15) from 0 to t to obtain 


n(Z,R,t) 


;k 



(A-16) 


where 

_2 

^ " 2(C{ + 2) 


and H denotes the step function 


H(t) = 


1 t > 0 
0 t < 0 


Equation (A-16) shows that at a given position Z all bubbles have the same 
radius and that bubbles arrive at height Z only after a time (Z/A)'/2. 

In particular, if h is the height of the thermal boundary layer then 
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$ 


n(h,R,t) 




«)'")"(.■(!)“). 


(A-17) 


This expression must now be used as a boundary condition for analyzing 
the bubble population in the second regime, i,e., above the boundary 
layer. 


A. 3. 2 The Bulk Fluid Above the Boundary Layer 

In the region above the boundary layer, the assumptions of the model 
are somewhat different. Here, the bubble no longer grows by evaporating 
more superheated liquid (since there is none), but rather by responding 
to changes in pressure and temperature as "it rises. Since the pressure 
is likely to vary as the bubble rises, it is no longer adequate to treat 
the bubble density pw as a constant, nor can it oe assumed that the bubble 
contains pure saturated vapor. In order to describe the bubble dynamics 
in this regime, a number of basic equations are needed. These are derived 
below. 
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The Mass Equation 


I 


For a fluid in general , we have 


3o 

3 t + V . pv 


0 . 


(A-18) 


(j* 


Integrate this equation over the bubble volume V(t), with surface 
S(t), to get 


J V • pv d^r = 0 . 

v(t) •'v(t) 

Use Leibnitz's and Gauss' theorems to transform 


(A-19) 


Uf pdV . f pVj . nj ♦ r 

\ V(t) •4(t) / Mt) 


pv • n^ d^S = 0. (A-20) 


Here, S(t) = surface, n^ = outer normal to surface, d^S * element of 
area on S(t), = velocity of surface point. 


Let 


M(t) = J pd^r = bubble 

V(t) 


mass . 


(A-21) 
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Then, 


dM 

Bt 


' r p(^ “ Jit) ‘ = evaporation rate. 

S(t) 


(A-22) 


The Energy Equation 
For a fluid in general 


Iy • P |u + ^ + gZ^ + V ' pv f h + --^ + gzj 


= V * (t • X " il) P9^2‘ 


(A-23) 


(Here, +Z = upward) Integrate this equation over V(t), then use 
Leibnitz's and Gauss' theorems to get 


p(u+4*9z)A-y p (. 

•'v(t) S(t) 


u + |-+gZjy^ 


+ /* P fh + |- + gzj V • n^ d^S = r 

4(t) 




• n^ d'^S 


S(t) 


•V(t) 


pv, dr. 


(A-24) 


Expect that kinetic energy terms (^) are small and that dissipation 
heat (/ I • • n^ d^S) is negligible. ‘^Then, since h = u + P/p, 
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g^/ pZd^r 


f pv^ dr - / 

'{t) ^(t) 


(A-25) 


Also, the potential energy terms cancel since 


PE.y 

•'M 


pgZ d r. 


Therefore, 


d'r = g^zff d^r pZy^ • n^ d^sj 


= g ^ J" ZV • py^ d^r + J" pZy^ • n^ 


- g V • py^Z d^r +J^ py^ • VZ d^r + pZy^ • n^ d^S^ 


y* pyZ • n^ d^S + pv^ d^r + pv^Z • n^ d^S^ 


9 pv^ d^r pZ (ys “ l) * JI5 d^S^ 


g pz d^r = g J* pv^ dV + g^* pZ (yj - y) • d^S. (A-26) 






Thus, the energy equation becomes 


■df pu d^r = / pu ( Vs - v) • n^ d^S - f 

•'v(t) 4(t) 4(t) 


Pv • n_ d S 
— — s 


/ 

-'s(t) 


J • n_ d S. 


(A-27) 


Take u and P uniform across V; let J denote last integral; then 


dt 


(Mu) = 


uj p (\ 


- v) • ru d^S 


-p/v-n, 


d‘"S - J. 


dM 


Since the first integral is ^ from tne mass equation. 


M 


du 

dt 




^ d‘"S - J. 


Under integral sign, take ±v^ to get 




^ d^S = / (vs - V) . n, d^S ^ ^ 


- / V, • a, 


d^S 


Since 


"/(t) 


^ dt 




d^S 


that is 



V 





(.Vs - v) 
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^ * 


then the energy equation becomes 

^ ^ = Pjf (Vj - »/) - Hg d^S - P ^ - j (A-28) 

The Momentum Equation 
For a fluid in general 

^ pv + V . pyv = - 7P + V . T - pgz . (2 = unit vector up) (A-29) 






Integrate over V(t) and convert integrals to get 


^ f PVdV -f ov 

J \l •'c 


V. • n. d S 
S S 


+ f P 11 ’ He 


d^S 


“ ■/ P d^S + / T . n^ d^S - gz / p d^ 


i.e. , 


^ r p V d^r = j 01 ( 1 ^ - 1 ) • Ug d^S -JPj}^ d^S 

•'v ‘'s *^5 

/ T . n 

= -s 


+ / T . d S - gzM. 


Note that 


/ 


P n d^S 

■~5 


buoyant force = - — M gz 


where is the liquid density and p the bubble density. Also, 
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t 



• d'^S = drag force = - 6TiuRyg 




(Here, refers to the mean bubble velocity.) 
i.e. , 


d 

dt 


J pv d^r = 
V 


/ 


pv (v^ - v) 


d^S ^ 


Mgz 


- eirpRVg - gMz. 


(A-30) 


In the first integral, approximate v^ by ys to get /ypv d^r = Mvg. In 
the second integral, replace the first y^ factor by yB- This is allowed 
if flow velocity into surface (due to condensation or evaporation) is 
small compared to mean bubble velocity. Then, 

ft J P (y^ - v) • d^S Mgz - S^iuRy^ - gMz. 

The integral is ^ from the mass equation so 


M - 67TuRy^ ) 
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Second Law of Thermodynamics 

In a system where M changes and heat dQ is transferred in, we have 


dS i ^ + sdM, i.e., dQ < TdS - TsdM. 


(A-32) 


Thus, for the bubble 

^ ■ dt - ' dt 3t 

or with S = sM, 

1 < T ds . dM 
' ^ ^ dt ■ =* dt 


-05tk||. 


Ideal Vapor Inside Bubble 

If specific heat is constant and s = 0 at T = T-j, P = P-| then 
P = Pv>"T 


\(PJ) = Cpyln (T/T^) - r In (P/P^) 


(A-35) 


= (Cpv - (T - T^). 


(A-36) 


This analysis also applies when there are small liquid droplets 
within the bubble. 


For reversible case, note J = 0 implies that s = constant. 

Pressure Discontinuity at Surface 

P (bubble) = P (liquid) + |2- (o = surface tension) (A-34) 
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Hydrostatic Pressure 


J 


P (liquid) = Pq " Pj^QZ 


(A-37) 


where 2 = distance above point where P = Pq. 


Geometry 


V=|ttR^. (A-38) 

Density 

M = pV. (A-39) 

1 

Displacement 


This is same as for boundary layer regime, i.e.. 



A. 4 SIMPLIFIED BULK FLUID MODEL 

Look at the simplest case: no mass or heat transfer = 0, J = 0), 
This does not rule out condensation of tiny droplets inside^the bubble. 

It just means that they stay in the bubble as mist rather than being 
absorbed into the surrounding liquid. 


( 
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We now have from mass equation [Equation (A-22)]: 









m 

dt 


0 


(A-41) 


with solution 

M (t) = Mq. (A-42) 

From second law ([Equation (A-33)] have ^ i 0* Approximate by 
reversible case ^ = 0, i.e.. 


s (t) = s (0). 


(A-43) 


Now, we can estimate the state of the bubble, which presumably leaves the 
boundary layer as saturated pure vapor. In rising into a lower pressure 
region at constant s, we must determine whether it superheats or, instead, 
changes to a mixture of saturated vapor plus droplets. In the former, 
pure vapor, case Equation (A-37) predicts that 


P T T 


'P/r 


2- 1- -fi- (1- f ) for T near T, 
r Iq I 


(A-44) 


If, instead, a vapor-liquid mixture results, then P is given by the satura- 

p 

tion pressure P^ (T). From corresponding state theory, ^ should be a 

^ T 

universal function of (c = critical point). In practice, this is 

I ^ 

approximate, and is given by the empirical fit*^ 


Pj (T) = P^ exp 


5.29 



(A-45) 


A-17 








Thus, 


where 


Then 


and 



(Tq. Pq at t = 0). 

f (T) = ^ 

T 

f(T) - T(Tq) + (T-Tq) f(TQ) + ... 



for T near Tq. 


(A-46) 


Thus, the coefficient of (l - is 5.31 | | for the liquid-vapor 

C V 0/ VqI 

mixture and p/r = 3/2 for the pure vapor (see Figure A-2). Cooling at 
constant s leads iimiediately into the liquid region, so the liquid-vapor 
case is the one that occurs. 



FIGURE A-2. Sketch of Pressure Temperature Data 
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Since hydrostatic pressure is assumed, we have [Equation (A-37)] 




Neglect surface tension so that 
P(Z) = Pq - Pl 92- 


Si nee the bubble is at saturation pressure, we have 


’’s(T) - >>0 - Pl 9 ; 
which determines T(Z), i.e.. 


" "^sat • 


If the corresponding states expression is used, this gives 


5.31 T 

" 5.29 - In/Pg - 


or with 


T = Tq at z = 0, 


5.29 + In 'c/p 
T(Z)lT„ \JJ\ 


i.29 + In /P_ V 


(A-47) 


Note that as L increases, T decreases. 


k 





To continue, the assumption of Ideal vapor gives 


= P/rT 


(Z) = P(Z)/rT(Z), 


Also, 


\ ■ S. ' 




For the liquid, assume Incompressibility so that 
Pl (Z) = Pl ■ const 

Sl (Z) ■ CpL In (T/^^) 
where T 2 is a reference temperature. 




Since 


With the above models we can construct the remaining properties: 


M|_ + = M 


Jz. + -L - \J 

Pl Py 


we have, after a little algebra. 


riM’t) 


(A-49) 
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Thus, 



I 

> t 


Now, since s is to be constant = Sq = Cp^ In ^^0^ r In (since the bubble 
is initially pure vapor), we must have ^ ^ 



which can be solved for p to get 


p(z) = r(z)-n - 


Pv 


(Z) Q (Z) 


Pvjzi\ 

^L(Z), 


SvVnp^ 


In In 

v^r 


in which 


- s, = Q/T 


with Q the latent heat of vaporization [q(Z) = Q [T(Z)]j. 


(if) 


(A-52) 
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Then 

where p^, = p,R at boundary layer interface. It remains to find the 

velocity v(Z) and the position 2(T). We note^^at this^^stage that u(Z) can 
be determined just as s(Z) was by taking u = jjp + jjp Thus, the 
energy equation should give no new information as can be seen if we 
rewrite Equation (A-28) using V = ^ and ^ ^ 9®^ 



Mn V du _ ^^0 V ^ 

“ dZ “ X dZ 

P 


(A-54) 


1 .e. , 


P d 


du 

az ■ 2 dZ 

P 




(A-55) 


which is satisfied identically since ^--^=T^=0. 

oz p dZ dZ 

To solve for v, rewrite Equation (A-31) using V = Mq/p and 


d 

at ■ 


V 


d 

dZ 


to get 


dZ N|j 



9 . 


(A-56) 


Since R(Z) and p(Z) are known. Equation (A-56) is of the form 
V ^ + f(Z) V = h(Z) 
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with f and h known. This is apparently not solvable in closed form,* 
f' even for simple f and h, although it can be solved numerically. For 

f. numerics, write as 

I = 

so that 

which is a nonlinear integral equation for v(Z). It can probably be solved 
successfully by iteration. 

Once v(Z) is determined, the position Z(t) follows from Equation (A-40): 


dZ 

Bt 


= v(Z) 


or 


or 


dZ 


dt 



dZ' 


t. 


With Z(t) determined, all other quantities available as functions of Z 
become explicit functions of t. 

It is disappointing that even this simplest case seems to require 
one numerical solution [v(Z)]; but perhaps it is surprising that all but 
one equation caji be integrated in closed form. 

With all quantities related to an individual bubble as determined 
above, the distribution function n(Z,R,t) can also be determined. Assuming 
r that no bubbles form by nucleation above the boundary layer, and th?w the 

! distribution of bubbles entering from the boundary layer at Z = h is given 

i- by Equation (A-17). 

i 

I 

i 


% 
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APPENDIX B 


DETAILED DEVELOPMENT OF MIXING MODEL 


This model is concerned with the high energy portion (above the 
boiling energy) of the critical mass curve. In this region the mixing 
mechanism, which is primarily high energy turbulence, is due to 
isotropic turbulence if 

L>>A 

where L = characteristic linear scale of the mixed system, and 

A = the size of the energy dissipating eddy. 

In this case the value of A can be related to power density from 
Kolmogoroff turbulence theory*^ 


A = 


3/4 

n 

172 




(B-1) 


where u = viscosity 
p = density 
P = power 
V = volume. 


A quantity, I., called the intensity of segregation can also be defined 
as^“ 


Is - e 


•t/y 


(B-2) 


where 


_ 

" ■ 6T 


The segregation intensity is also the fraction of one component actually 
not yet mixed into the other component. 

Define a quantity "degree of mixing" 

f = ^ 

^ M 

where m = amount of A mixed with B, and 
M = total amount of A in the system. 

Then f = 1 - 1^ 
or 


i 


M 


= 1 


-t/y 


(B-3) 


Equation (B-3) describes the mixedness as a function of time for a batch 
process. 

B-1 




Now consider a turbulent mixing system in which there is a con- 
tinuous influx of component A given by dM/dt. Each increment can be 
thought of as undergoing batch mixing over a time period t. That Is 


dm » dM ^1 - e^^ ‘ 


(B-4) 


where dM is added at time t'. Equation (B-4) must be integrated over 
all dm to get total mixing as function of time. 




(B-5) 


In the special case where M is dropped through a hole from a 
constant height (the ullage space) it is reasonable to view the flow 
as constant velocity since the head and opening resistance are approxi 
mately constant. Therefore, 


dM = k^df 

wnere is a constant. Since m(o) = 0 

M = K^t'. 

Let <2 = 1/y. then Equation (B-5) becomes 


m = /* 1 - K^df 

*'o L 


after substitution of Equation (B-6). After integration 


Expand (l - e'*^2^) in a series: 


2 2 3 3 4 A 

„ . Vt 

^ -1/4. A J. ^ C 

1 - e 2 - < 2 t - + 3 ^ - -y| 


(B-6) 


(B-7) 


(B-8) 


(B-9) 
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then 


m 




for small t 


m 


K, K,t' 
% I 2 


K^K2^t^ 

-^T — “Tl — 

• • 

Kjt^ 


(B-10) 


Now consider the power input. The system gains energy from the 
falling fluid and, if there are no appreciable friction losses, the 
energy is 


E = Mgx 

where x = ullage space height. 
Power is 


dE _ dM /„ \ 

dt - dt (9V 

since g,x are constants. 

Since ‘ K-j = a constant 

then ^ = K^gx = a constant 

and the coefficient K 2 is a constant given by 

6p^K.jgx 

■^2 % 5 / 2 ^ 1/2 • 

If 

E = Mgx = K^tgx = K^t 


(B-11) 


(B-12) 
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then 


or 




(B-13) 

♦ ■ ( 


The validity of the approximation in Equation (B-10) was tested 
numerically. Typical data in the Farber studies’ gives mixing times 
of about 1 sec with yields corresponding to 3000 lb in a 25,000-lb 
system. If 25,000 lb are added in 1 sec, then = 25,000. K 2 can then 
be estimated from Equation (B-10) by trial and error; in this case 
K 2 = 0.26. Using Ki , K 2 as calculated. Table B-1 was generated. For 
these values of K] and K 2 (which are of the order of magnitude expected) 
the approximation using one term is very good up to 1 or 2 sec. At 
most, the approximation is conservative in that larger mixed masses are 
predicted than are calculated by the exact equation. 


TABLE B-1 . Comparison of Exact Expression 
[Equation (B-9)] with Series 
Approximation 


Time, 


The Mass Mixed 

(M), lb 


sec 

Exact 

1 Term 

2 Terms 

3 Terms 

0.1 

32.2 

32.5 

32.2 

32.2 

0.2 

127.8 

130.0 

127.7 

127.8 

0.3 

2P5.0 

292.5 

284.9 

285.0 

0.4 

502.4 

520.0 

502.0 

502.4 

0.5 

778.4 

812.5 

777.3 

778.4 

0.6 

1111.5 

1170.0 

1109.2 

1111.5 

0.7 

1500.1 

1592.5 

1495.9 

1500.3 

0.8 

1943.0 

2080.0 

1935.8 

1943.3 

0.9 

2438.6 

2632.5 

2427.2 

2439.2 

1.0 

2985.7 

3250.0 

2968.3 

2986.6 

2 

11011.6 

13000.r 

10746.7 

11039.6 

5 

55051.1 

8125U.0 

46041.7 

57*34.4 
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QUENCHING DISTANCE IN DIELECTRICS 


Quenching distances discussed in this report are values measured by 
experiment using metal electrodes. The values are a measure of the 
greatest wall distance for which wall dissipation of energy or ions can 
inhibit combustion. A quenching distance so measured is approximately 
the electrode spacing for which spark energy needed for ignition is 
smallest. For electrode spacing closer than the quenching distance, 
wall scavenging of ions and energy makes a more energetic pulse necessary 
to start combustion. 

This study considers ignition by discharges in a dielectric environ- 
ment, e.g., between dielectric liquid droplets or across bubbles in 
dielectric liquid. The dielectric should be less effective in suppressing 
ignition than the metal electrodes for two reasons: 

1) Ions are less strongly attracted to the dielectric wall than to 
metal . 




2) The lower thermal conductivity of the dielectric wall results in 
slower conductive heat dissipation of energy deposited on the surface 
by radiation, and more re-radiation into the reaction region. 

The ion dissipation with low wall distances maiy occur by the follow- 
ing sequence of events. Following a discharge, ions and electrons are 
both accelerated toward walls by electrical image forces. Electrons 
arrive. in significant numbers first because of their more rapid accelera- 
tion. Following a loss of space charge neutrality, however, ions deposit 
rapidly because of the stronger collective field force (for short distances 
to walls). The quenching distance is on the order of the distance cleared 
of reacting ions in the mean ion reaction time in the combustion process. 


A particle of charge q a distance s from a metal wall experiences an 
image force from induced surface charge given (in rationalized MKS units) 
by 


F M = _ _J gi_ 

S 4ireo (2s)^ 


(C-1) 



i 


1 


I 




y 


The corresponding force f om induceo polarization of a dielectric wall 
distance s from the particle is 


P D ^ _ _1_ M 
s ■ “ Atveo (2s)^ 


(C-2) 


K 

k 





where K is wall dielectric constant. The negative sign in Equations (C-1) 
and (C-2) comes from measuring distance s positive from wall while the 
force is directed toward the wall. Neglecting colli sional damping and 
forces from other particles, the particle will be accelerated toward 
the wall according to 


m 


d2s 



(C-3) 


where c is a product of factors in E-1 or E-2. 


Multiplying Equation E-3 by ds/dt gives 

d^s ds _ c ds 

dt2 dt " ■ ^ 3t 



(C-4) 


Integrating Equation (C-4) from time t=o at which ds/dt is zero and s is 
Sp to a later time t gives 


1 

2 




(C-5) 


Taking the square root (and noting ds/dt is negative) and rearranging 
gives 



and integrate Equation (C-6) from position s at time t=o to position 
s=o, i.e., hitting wall at time t=T: 


C-2 


' ' . 






(^0^^ - ^0^^) 
-s^ 


1/2 


1 1 / 2s„s - s„ 

♦ ,2 '— Sin-’ ' ° “ 


s=o 


s=s, 


iff ‘ 


t=T 

t=0 


(C-8) 


or 


1 c 3/2 _ /2cy/2 . 

2 ^0 - "Ui '• 


(C-9) 


Isolating initial distance Sq in Equation (C-9) gives 


=0 = iifeif 


(C-10) 


The quenching distance should be on the order of the distance from the 
wall cleared of ions in a mean ion reaction time. Setting T to the mean 
ion reaction time on the right hand side of Equation (C-10) should give an 
approximate quenching distance on the right hand side. We assume that the 
migration time to the wall of the faster charge species sets the range of 
quenching. The value of c is different for the metal and the dielectric. 
Hence from Equation (C-10) 



For dielectric constants K in the range 1.25 to 1.5, the quenching 
distance by (C-11) should be 0.5 to 0.6 of those for metal. 

The foregoing development made no mention of Debye-Huckel screening 
in the interaction of charges with walls. At the threshold for ignition, 
the free electron number density should be quite low, and hence the screen 
ing distance should be long. The thermal ignition temperature for 
hydrogen-oxygen mixtures is found to be about 585°C. The Saha equation's 
predicts a low level of thermal ionization at this temperature. 
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APPENDIX D 

CONTACT POTENTIAL FUNDAMENTALS 


Each homogeneous object has Its characteristic Fermi level which is 
a function of temperature. The Fermi energy level is the electron energy 
level above which the probability of occupation of an electron state is 
less than 50% and below which the probability of occupation is greater 
than 50%. Specifically, the probability of occupation of an electron 
state of energy e is 


f(e. T) = L zr (D-1) 

1 + exp |^(e-M)/kTj 

where y is the "chemical potential" of an electron or the Fermi energy 
level. The reference point of energy must be the same for both e and y. 
The form for f says that electron states of energy e much less than y 
for- which 


.1 


I 


^ « -1 (D-2) 

are filled with a probability only slightly less than one. Electron 
states of energy greater than y for which 

>> 1 (0-3) 

will have a probability much less than one of being occupied. 

When two dissimilar objects of homogeneous material are brought 
into contact, electrons tend to flow so that states of equal energy 
(relative to the same reference point of energy) have the same proba- 
bility of being filled in both objects. To see this, consider two 
metals at temperature T whose electron energy state occupancy is shown 
in Figure D-1. The electrons in Material 1 see a potential well tending 
to bind them within the metal. The most tightly bound electrons of the 
conduction band have an energy egi relative to the energy of a distant 
isolated electron in vacuum. Typically, ebI (binding energy) will be a 
negative energy of the order of 10 eV. The conduction electrons in 
Metal 1 will occupy states which have the energy eg] plus varying amounts 
of kinetic energy up to a kinetic energy of approximately e|(p-| (kinetic 
energy at the Fermi level of Material 1). The transition from occupied 
to unoccupied states occurs in a region of width -vkT about this value. 

The Fermi level for isolated Material 1 we give relative to that of an 
isola d distance electron at rest in vacuum as 




D-1 



I 


^ ^B1 ^KFl. 


(D-4) 


Note that -vii is the energy needed to remove the most energetic electrons 
from highly occupied states into vacuum from the neutral metal; note also 
that our reference point for Fermi energy may differ from that used 
elsewhere. 


EMPTY 

STATES 

PARTIALLY 

FILLED 

FILLED 

STATES 



METAL 1 • 


ENERGY OF AN ISOLATED 
ELECTRON AT REST 
IN VACUUM 

FERMI LEVELS 
(505fc OCCUPANCY) 


ENERGIES OF MOST TIGHTLY 
BOUND ELECTRONS 
IN CONDUCTION BAND ^ 



METAL 2 


FIGURE D-1 . Energy Level Schematic for Two Isolated 

Neutral Pure Metals, (Note that we measure 
energies positive upward, so egi , eB2. ui , 
and u2 negative numbers.) 


When the metals from Figure D-1 are brought into contact, electrons 
will initially flow from Metal 1 to Metal 2 (because ui > ^ 2 )- 
infinitesimal fraction of the electrons need to flow, however, to create 
an electric field which stops the flow. The well depths of the two 
metals will have been shifted by the electrostatic contributions qeV] 
and qeV 2 » respectively, when equilibrium is established. Here, qe is 
the charge on an electron. The equilibrium situation for the bulk mate- 
rials is shown in Figure D-2. The condition for no net electron flow 
between the two media is 


y? + qg Vi = + q^ (D-5) 

Here, the superscript zero on m? and merely indicates that it does 
not include the electrostatic potential contribution which we have 
included in our definition of u. For an inhomogeneous medium at uniform 
temperature we are led to postulate an electric current density 

0 

J = -0 V (^ + V) (D-6) 


D-2 





METAL 1 


METAL 2 


FIGURE D-2 . Energy Level Schematic for Two Pieces of Dissimilar 
Metals Brought into Contact 

in the absence of magnetic field. This reduces to the familiar J = -a VV or 
J = oE with conductivity a for a homogeneous medium. This is in agreement 
with Landau and Lifshitz.^° With nonuniform temperature, a current contri- 
bution proportional to temperature gradient must be added: 

/ 0 \ 

j = o I- V(^ + V) - aVTl. (D-7) 

For a contact between dissimilar semiconductors or good resistors, 
similar considerations apply. A charge transfer occurs on contact to 
equalize the Fermi levels in the two materials, although Fermi levels 
will for resistors or semiconductors typically be mathematical energies 
only [in the Fermi-Dirac distribution. Equation (D-1 )] , not energies of 
actual allowed electron states. Also the approach to equilibrium will 
be very slow compared to that in meta’.s. 

The potential difference V 2 - V] between the dissimilar materials 
comes phyically from a charge distribution at the interface. This 
interface charge distribution for good conductors is very thin, and is 
frequently considered as a dipole layer. For a germanium junction with 
10^5 donor atoms per cubic centimeter in contact with a metal, the space 
charge region is estimated to be lO"^ cm thick^^ at room temperature. 

For electrolytes as well as semiconductors, the thickness of the space 
charge region varies with the inverse square root of the number density 


♦ *w 

« P' 

; 

; f 

/ 




of current carriers. This finite thickness of the space charge region 
(particularly at cryogenic temperatures) for low temperature materials 
with an extremely low density of charge carriers gives an opportunity 
for pronounced charging effects in vaporization Induced bubbling and 
frothing. It should be noted that the slow relaxation time in highly 
resistive liquids or vapors implies that the charge separation may 
occur only if the vapor bubbles have traveled through the medium long 
enough for charge separation to have occurred. 
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